Structure effect on one neutron removal reaction using relativistic mean field densities in Glauber 

model 
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We calculate the one neutron removal reaction cross-section (cr_i n ) for few stable and neutron-rich halo 
nuclei with 12 C as target, using relativistic mean field (RMF) densities, in the frame work of Glauber model. 
The results are compared with the experimental data. Study of the stable nuclei with the deformed densities 
have shown a good agreement with the data, however, it differs significantly for the halo nuclei. We observe that 
while estimating the <J-i n value from the difference of reaction cross-section of two neighboring nuclei with 
mass number A and that of A-l in an isotopic chain, we get good agreement with the known experimental data 
for the halo cases. 
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I. INTRODUCTION 

From last two decades or more, the exploration of neutron- 
rich nuclei is an important branch in Nuclear Physics research. 
It is a source of observance of new phenomena and dynam- 
ics. This is possible due to the development of accelerator 
techniques for Radioactive Ion Beams (RIBs) in various lab- 
oratories around the globe. Experimental methods and theo- 
retical analysis have been widely used to collect information 
about the structure information, such as nuclear size, valence 
nucleon distribution and halo structure of these exotic nuclei. 
The measurement of various reaction observables like total re- 
action cross-section a r , one- and two- nucleon removal cross- 
section (a -i n , (T-2n) and the longitudinal momentum distri- 
bution P| | are some of the established quantities for such stud- 
ies. 

The relativistic mean field (RMF) or the effective field the- 
ory motivated RMF, i.e., the E-RMF provides the internal 
structure or sub- structure information of the nuclei through 
the density distributions, using as an input while calculating 
the observables in conjunction with Glauber model lll|2|]. A 
systematic study of one- and two-neutron knockout data for 
15-i9q explained beautifully while using shell model, which 
gives a consistent structure information not only for the sta- 
ble nuclei but also neutrons at the boundary [3]. It is well 
known that narrow fragment momentum distribution reflects 
large space distribution of the valence nucleon and there is a 
correlation of the magnitude of the a_i n with the width of 
the P\ | , in approaching the nucleon (neutron or proton) drip- 
lines. However, one-neutron removal reaction cross-sections 
provide important nuclear structure information complemen- 
tary to that obtained from P\ | . 

In the present paper, our aim is to calculate the a r and a_i n 
by using densities obtained from the RMF and E-RMF for- 
malisms 111 |2|] in conjunction with Glauber model. It will be 
shown for the nuclei near the drip-line which possess a halo- 
structure fails to be explained by the standard evaluation of 
a_i n . Contrary to this estimation, the difference in total re- 
action cross-section between two consecutive neighboring nu- 
clei in an isotopic chain better matches with the experimental 
data. 



II. THEORETICAL FRAMEWORK 

The use of RMF and E-RMF formalisms for finite nuclei as 
well as infinite nuclear matter are well documented and details 
can be found in |4j,[5|] and (aLZD respectively. The working ex- 
pressions for density profiles and other related quantities are 
available in 01 SB E 0]- The details to calculate cr r us- 
ing Glauber approach has been given by R. J. Glauber 0]. 
This model is based on the independent, individual nucleon- 
nucleon (AW) collisions in the overlapping zone of the col- 
liding nuclei, and has been used extensively to explain the ob- 
served total nuclear reaction cross-sections for various sys- 
tems at high energies. The standard Glauber form for thetotal 
reaction cross-sections at high energies is expressed as | 



cr r = 2tt J b[l - T(b)]db , 



(1) 



where T(b) is the transparency function with impact param- 
eter b. The function T(b) is calculated in the overlap region 
between the projectile and the target assuming the interaction 
is formed from a single NN collision. It is given by 



T(b) = exp 



• a ij / ds Pu ( s ) Ppj ( 



• (2) 



The summation indices i and j run over proton and nucleon 
and subscripts p and t referred to projectile and target, re- 
spectively. The experimental nucleon-nucleon reaction cross- 
section aij varies with energy. The z-integrated densities 
p(uj) are defined as 



p(u) 



dz 



(3) 



with u 2 = x 2 +y 2 . The argument of T(b) in Eq. (2) is — s 

which stands for the impact parameter between the i th and j th 
nucleons. 



The original Glauber model was designed for high en- 
ergy approximation. However, it was found to work reason- 
ably well for both the nucleus-nucleus reaction and the dif- 
ferential elastic scattering cross-sections over a broad energy 
range 1 10]. To include the low energy effects of NN interac- 
tion, the Glauber model is modified to take care of the finite 
range effects in the profile function and Coulomb modified 
trajectories lUlltil . The modified T(b) is given by (HQ, 



T(b) exp 



" f p [J2 [ r « * + *)] P P i (*) P tj (S) dsdi 



The profile function (6 e //) is defined as HI] 



(4) 
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FIG. 1: The spherical proton (p p ) and neutron (p n ) density obtained 
(5) from RMF (NL3) parameter set for various isotopes of (a) Carbon 
and (b) Boron. 



with b e ff = b — s + t , s and t are the dummy variables 
for integration over the z-integrated target and projectile den- 
sities. The parameters a at at, a at at, and @nn are usually 
case-dependent (proton-proton, enutron-neutron or proton- 
neutron), but we have used the appropriate average values 
from Refs. i9l fl . The deformed or spherical nuclear den- 
sities obtained from the RMF and E-RMF models are fitted 
to a sum of two Gaussian functions with suitable co-efficients 
Ci and ranges en chosen for the respective nuclei which is ex- 
pressed as 



p(r) = ^Cjexpl-ajr 2 }. 



(6) 



Then, the Glauber model is used to calculate the total reaction 
cross-section for both the stable and unstable nuclei consid- 
ered in the present study. 

The expression for one nucleon removal reaction cross- 
section <7_i n (7) is given by II 1211 



E 



d k a, 



a=(/c,g=0),C5 



(7) 



where o- a =(k,g=Q),c are the possible final states ac. In the 
present formalism, it is considered that the projectile nucleus 
breaks up into a core and the removed nucleon. The core C 
has an internal wave function <p g and the one-nucleon, i.e., the 
departed nucleon has an asymptotic momentum fik in the con- 
tinuum state with respect to the core. The core is considered 
to be in the ground state (g = 0) at the time of the collision. 
The total cr_ ln (J) can be separated to an elastic (cr£? ln ) with 
c = and inelastic (cr!f( n ) part having c as non-zero. The 
<j e \ n and (J l l\ n is expressed as [QT 



<r-m(I) = J dh i< ^ I e 



- < 0o I e 



2Imxct(b c )-2Imx-int(b c +s) | 
iXct{b c )+iX-int{b c +s) | ^ > |2j 



J db{< 0o 



-2Imxct(b c ) 



-2Imxct{bc)-2Imx-int{b c +s) 



k >},(9) 



here Xpt is the phase shift function and 0o is the valence wave 
function (the wavef unction of the removed nucleon). The 
notation and the numerical procedure of calculation of one- 
nucleon removal reaction cross-section are followed from Ref . 



III. CALCULATIONS AND RESULTS 

We obtain the field equations for nucleons and mesons from 
the RMF and E-RMF lagrangian. For deformed case (RMF 
only), these equations are solved by expanding the upper and 
lower components of the Dirac spinners and the boson fields 
in an axially deformed harmonic oscillator basis. The set of 
coupled equations are solved numerically by a self-consistent 
iteration method taking different inputs of the initial defor- 
mation /?o iBIlH]. For spherical densities both for RMF and 
E-RMF models, we follow the numerical procedure of Refs. 
0]. In our calculation, the constant gap BCS pairing is used 
to add the pairing effects for open shell nuclei. The centre- 
of-mass motion (cm.) energy correction is estimated by the 
usual harmonic oscillator formula E c . m . = |(41t4 -1 / 3 ). 

Comparing the binding energy (BE) of the calculated solu- 
tions, the maximum BE and the corresponding densities [p p 
(proton) and p n (neutron)] are for the ground state. All other 
solutions are the excited intrinsic state including the spherical 
one. Since the main input in the Glauber model estimation is 
the RMF or E-RMF densities, it is important to have an in- 
formation of these quantities. We have plotted the spherical 
>p p and p n for both proton and neutron of Carbon and Boron 
isotopes in Figure 1 using RMF (NL3) parameter set Q3|] . As 
(Expected, we find an extended density distribution for proton 
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compared to neutron in case of 9 C and 8 B due to the proton- 
rich nature of these two nuclei. The value of p n and p p are 
almost similar for 12 C which can be seen from Figure 1. Ex- 
tension of p n is much more than p p for rest of the nuclei. It 
is maximum for 19 C and 15 B in Carbon and Boron isotopic 
chains, respectively, because of the high neutron to proton ra- 
tio for these cases. 

In the present study of a r and cr_i n , first we use the spher- 
ical density obtained from RMF (NL3)|H and E-RMF (G2) 
B]. The results are presented in Table 1 for 9,12,13,15,17,19^ 
and M2, 13,14,15^ isotopes with 12 C— target at various pro- 
jectile energies. These results deviate considerably from the 
data IT 61-1 1911 which are quoted in the table. For example, in 
case of 9 C+ 12 C, the observed value of a_i n is 48 ± 8 mb as 
compared to the estimated results of 8 1 and 96 mb with NL3 
and G2 parametrization, respectively. Note that the <j\ n for 
8 B -f 12 C and 9 C -f 1 2C systems are one-proton removal re- 
action cross-section, which may be followed throughout the 
text and Tables. However, in rest of the systems, a\ n will be 
refered as one-neutron removal reaction cross-section. Sim- 
ilar discrepancy is also seen for other cases. A further in- 
spection of the table shows that the experimental one-neutron 
removal reaction cross-section for some selected cases coin- 
cide well with the prediction. We also used the method of B. 
Abu-Ibrahim et al. 1 12] to calculate the one-neutron removal 
reaction cross-section cr_i n (/7), which obtain by the differ- 
ence of total reaction cross-section of two neighboring nuclei 
with mass number A and A-l in an isotopic chain. This pre- 
scription is suitable only for halo projectile and may not be 
applicable for general cases. The expression is given by |H2II : 

<7_i n (/I) = a r ( A Z) - a r { A - x Z), (10) 

and the values are inserted in Table 1 for comparison. The 
a-i n (II) differs significantly from the experimental data for 
all the cases. It is important to recall that the effect of deforma- 
tion is nominal in the evaluation of <j r which is reported in our 
earlier publications 0112]. In these papers, the Glauber model 
with RMF (NL3, NL-SH) and E-RMF (G2) densities show 
a good agreement with experimental data for a r and elastic 
differential scattering cross-sections da/dQ, which in general 
justify the model independency of the calculation with various 
relativistic parametrizations. 

Unlike to the total reaction cross-section, the a_i n obtained 
from the Glauber model, depends very much on the structure 
information of the projectile and target nuclei, i.e., input den- 
sities of these systems. 

The nuclear single-particle energy e nyP for the last occu- 
pied orbit is very important for a reaction process. Thus, it 
is worthwhile to analyse the e n?p of the valence nucleon of 
the projectile and target nuclei. For simplicity, the spheri- 
cal single-particle energy for the last occupied orbit for pro- 
ton e p and neutron e n with RMF (NL3) and E-RMF (G2) are 
compared. As expected, a small variant in e p or e n makes a 
remarkable change in cr_ ln (JJ) [Eq. (10)] for many cases. 
For example, the one neutron removal reaction cross-section 
a- ln (II), for 12 C -f 12 C are 55 and 39 mb for RMF (NL3) 
and E-RMF (G2) with their single-particle energies e p = 
-15.66, e n = -18.97 and e p = -13.16, e n = -16.16 MeV, 



respectively. This discrepancy is minimum in the calculation 
of <7_i n (J) [Eq. (7)] with a lone exception for 13 B -f 12 C 
system. As the valence e n ^ p plays a major role to determine 
the reaction observables, one needs to reproduce these values 
with the experimental observation. This can be achieved by 
a small adjustment of the parameters in the relativistic mean 
field formalisms. However, the philosophy of RMF or E-RMF 
of single set of parametrization for the entire domain of nu- 
clear landscape goes against this parameter fiddling. Keeping 
this in mind, the quality of the results is compromised slightly 
using with the original values of NL3, NL-SH or G2 sets. 

Apart from the single-particle energy, the structure effect of 
the participating nuclei is crucial for a reaction study. In this 
context, it is interesting enough to see this effect (deformation 
effect) on cr_i n . We repeat the calculations for cr_ ln (J) and 
a-i n (II) with the deformed densities (RMF only) as input in 
the Glauber model [ 1]. We obtain spherical equivalent of the 
axially deformed densities using equations (3) and (6) follow- 
ing the prescription of Refs. (ll0]. The NL-SH parameter set 
[ 20] for this purpose is used and the results are listed in Table 
2. The reason to change the NL3 to NL-SH is the unavail- 
ability of converged ground state deformed solution with NL3 
for very light mass nuclei, which is a situation in the present 
study 1 21]. Also, the NL-SH parametrization is reasonably a 
better parameter set and we expect similar outcome from this 
force parameter. 

Due to similar reason as mentioned for the spherical nu- 
clei, the deformed densities for some selected cases are im- 
perative to analyse. Our earlier work on density study sup- 
ply us enough signature about the complicated sub- structure 
Il22ll . The clustering and sub-structure of these deformed neu- 
tron and proton density distributions are demonstrated in Fig- 
ure 2. The density contours presented are in boxes of width 
and height 6 fm. A uniform contour spacing of 0.01 fm -3 
is used for proton and neutron densities. The z-axis is cho- 
sen as the symmetry axis, the densities are evaluated in the 
zp plane, where x = y = p. In ref. 112211. it is noticed that 
12 C possesses a 3a— cluster with a tetrahedral configuration. 
The same structure is reproduced in the present study with an 
oblate shape. The structure of the neutron deficient 9 C nu- 
cleus has a prolate ground state and that of the neutron-rich 
19 C has an oblate ground state deformation. In all the three 
cases, the density plot informs that the central part of the nu- 
cleus is a compact core, which is surrounded by a thin layer 
of nucleons. The structure of the internal core for both pro- 
ton and neutron have different density distribution from 9 C to 
19 C. The shape of 19 C proton density distribution looks like 
a perfect dumb-bell. Thus, it has a maximum probability to 
to have the structural effects on neutron removal reaction. On 
the other hand, the total nuclear reaction cross-section is less 
influenced by deformation, may be because of the averaging 
in input density in the Glauber model calculations. 

The results obtained from the deformed densities are tabu- 
lated (see Table 2). Table 2 shows that most of the cr_i n (7) 
[obtained from Eq (7) 1 matches quite well with the exper- 
imental data of Ill6l4l9fl and few of them do not agree. On 
the other hand a_i n evaluated from Eq. (10) coincide with 
only 19 C -f 12 C experimental data lH3l . Of particular interest 
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TABLE I: One-neutron- for 12 > 13 > 15 > 17 > 19 c an d one-proton removal reaction cross-sections cr_i n (in mb) for 9 C and 8 B with 12 C— target 
using spherical RMF and E-RMF densities obtained from NL3 and G2 parameter sets, respectively. The available experimental data fllq - [l9n 
is given for the comparison. The last occupied proton and neutron single-particle energies are also given. 



Projectile 


Energy 


CT-ln 


G- 




a-m(II) 


e P (MeV) 


e n (MeV) 






Exp. 


RMF 


E-RMF 


RMF 


E-RMF 


RMF 


E-RMF 


RMF 


E-RMF 




285 


48(8) 


81 


96 


34 


31 


-4.07 


-4.15 


-13.81 


-13.90 


12 r i 

c 


1050 


44.7(3) 




37 


55 


39 


-15.66 


-13.16 


-18.97 


-16.16 


13 C 


800 




29 


28 


47 


37 


-16.81 


-15.30 


-17.77 


-16.35 


15 C 


54 


137(16) 


95 


88 


37 


38 


-19.18 


-20.47 


-10.57 


-10.39 


17 c 


904 


129(22) 


101 


95 


45 


32 


-23.93 


-22.65 


-11.53 


-11.57 


19 C 


910 


231(51) 


134 


128 


35 


37 


-27.25 


-25.89 


-12.61 


-12.50 


8 B 


285 


89(2) 


103 


107 


47 


33 


-23.16 


-20.73 


-31.50 


-28.96 


12 B 


67 


81(5) 


54 


49 


47 


44 


-39.95 


-34.18 


-14.20 


-13.49 


13 B 


57 


59(4) 


54 


38 


46 


43 


-39.88 


-35.49 


-7.81 


-7.85 


1A B 


50 


153(15) 


70 


64 


48 


38 


-41.00 


-36.85 


-8.40 


-8.41 


15 B 


43 


108(13) 


87 


89 


51 


39 


-42.16 


-38.23 


-8.77 


-8.96 



TABLE II: Same as Table 1, but using both spherical and deformed NL-SH densities for 9 ' 1215 ' 1719 C and 8 ' 12131415 J B projectiles taking 12 C 
as target. The available experimental data is also displayed for the comparison lfT^ - fl9h . The last occupied neutron and proton single-particle 
orbits are also given and it is denoted by the Nilsson index [NnsAjQ 71 . 



Projectile 


Target 


Energy 


(J-ln 


cr-in(I) 


CJ-ln 


(II) 


ft 


neutron 


proton 








Exp. 


Sph. 


Def. 


Sph. 


Def. 




[Nn 3 A] 


e n (MeV) 


[Nn 3 A] 


e p (MeV) 


9 C7 


V2 C 


285 


48(8) 


70 


66 


24 


03 


0.36 


[110]1 - 


-16.122 


[101]3 _ 


-2.772 




12 C 


1050 


44.7(3) 


39 


49 


45 


293 


-0.21 


[101]1 - 


-16.952 


[ioi]i- 


-13.752 


15 C 


12 C 


54 


137(16) 


92 


130 


35 


15 


0.25 


[220] 1 + 


-2.611 


[101]3 _ 


-19.631 


17 c 


12 C 


904 


129(22) 


122 


120 


31 


32 


0.45 


[211]3 + 


-2.705 


[101]3 _ 


-21.906 


19 c 


12 C 


910 


233(51) 


134 


152 


31 


263 


-0.43 


[202] 3 + 


-3.316 


[ioi]i- 


-27.048 


8 B 


12 C 


285 


89(2) 


101 


99 


34 


10 


0.63 


[lioji - 


-14.054 


[101]3 _ 


-1.880 


12 B 


12 C 


67 


81(5) 


41 


65 


45 


70 


0.18 


[101]1 - 


-6.212 


[101]3 _ 


-15.585 


13 B 


12 C 


57 


59(4) 


53 


46 


44 


16 


0.10 


[noji- 


-7.378 


[101]3 _ 


-18.321 


14 B 


12 C 


50 


153(15) 


67 


97 


32 


15 


0.38 


[220] 1 + 


-1.992 


[101]3 _ 


-18.686 


15 B 


12 C 


43 


108(13) 


88 


104 


68 


73 


0.59 


[220] 1 + 


-2.611 


[101]3 _ 


-19.631 



amongst the nuclei investigated here are B and C which 
based on the relatively weak binding of the valence neutrons 
and measurements of one neutron removal cross-sections have 
suggested to be one neutron halo systems 12311 . The single- 
particle energy for proton e p and neutron e n for the last oc- 
cupied orbit are given in the 11th and 13th column of Table 
2. The last proton for 8 B and 9 C and the outer most neutron 
for 14 ' 15 £? and 17 ' 19 C are loosely bound which are the pos- 
sible candidates for either proton-halo (or skin) or neutron- 
halo (skin). Going back to the analysis of Figure 2, which 
gives us enough indication for the absence of halo-like struc- 
ture in 9 ' 12 C. Contrary to the case of 9 C and 12 C a thin-layer 
of neutron distribution spread spatially to a large extent in 
case of 19 C, which looks like a halo-nucleus. This behav- 
ior is also reflected in the one neutron removal reaction cross- 
section. In this particular case of 19 C, the calculated result 
a-i n (II) = 263 mb is more closure to the experimental value 
of a -i n = 233 ±51 mb than the a_i n obtained by using Eqn. 
(7). If we recall the statement of Abu-Ibrahim et al. 1 12], then 
19 C is a halo-nucleus. On the other hand, the measured one 
neutron removal reaction cross-section of 14 B is larger than 
its neighbors suggest the weak binding of the last neutron and 



extended valence density distribution. 

Summarising the whole discussions of Tables 1 and 2, in 
general, one can say it very briefly that except few cases like 
12,17 C+ 12 C and 8 ' 13 B+ 12 C the spherical density used from 
RMF (NL3) and E-RMF (G2) fails to reproduce the data. 
When we use the deformed densities to evaluate the one neu- 
tron cross-section, the predicted cr_ ln (J) matches reasonably 
well with the experimental measurement. In this case, only the 
result of the system 19 C+ 12 C deviate from the observation. 
However, the results predicted by Eq. (10) disagree largely 
with the experiments irrespective of the densities used. Com- 
plementating to Eq. (7), the <r ln (7T) matches with the lone 
case 19 C+ 12 C agreeing with the prediction of Ref. II 1 211 . 

In Figure 3, we have presented the cr_i n (7) with various 
incident energies for 19 C + 12 C using the spherical NL3, NL- 
SH and G2 densities in the Glauber model calculation. We 
also compare our results with the deformed NL-SH densi- 
ties obtained from the axially deformed RMF. All the spher- 
ical densities reproduce similar elastic (cr^ ln ) and inelastic 
(0"-in) one neutron removal reaction cross-section. The de- 
formed NL-SH densities has a large impact on the evaluation 
of a -i n unlike to the total nuclear reaction cross-section cr r , 
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P = -0.283 





p =-0.288 



19 C 




FIG. 2: The axially deformed density distribution for 91219 c with 
RMF (NL-SH) parameter set. 
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FIG. 3: The energy dependence of the neutron removal cross-section 
for 19 C+ 12 C system using spherical densities of RMF (NL3), RMF 
(NL-SH) and E-RMF (G2) parameter sets for both elastic and in- 
elastic processes. The result obtained by deformed RMF (NL-SH) 
densities is also given for the comparison. 



which is evident from Figure 3 and consistent with cr_ ln (J). 
The deformed (J % l\ n value constantly over-estimate the spher- 
ical cr l ^[ n starting from low to very high incident energy of the 
projectile. 



IV. SUMMARY AND CONCLUSION 

In summary, one neutron removal reaction cross-sections 
for the neutron-rich isotopes have been calculated using the 
densities obtained from RMF (NL3) and E-RMF (G2) for- 
malisms for spherical and deformed NL-SH parameter sets. 
The dependence of a_i n on single-particle energy of the last 
occupied nucleon is seen in our present calculations. That 
means, although the total nuclear reaction cross-section does 
not show a significant difference, the a_i n values differ from 
each other depending on the NL3 or G2 parameter set. The 
a -i n are in good agreement with the experiments, when we 
consider the deformation effect in the densities. The Glauber 
model fails for halo systems and in this case 19 C+ 12 C is a typ- 
ical example. In such case, the difference between the total re- 
action cross-section from the consecutive nuclei is applicable 
to evaluate a_i n (II). It is also concluded in the present paper 
that the deformation effects for one neutron removal cross- 
section is very much crucial unlike to the total reaction cross- 
section a r . In other words, the Glauber model reproduce the 
experimental data pretty well while considering the deformed 
densities for stable nuclei as projectile. On the other hand, 
when we estimate the difference of reaction cross-section of 
nuclei with mass number A and that of A-l in an isotopic 
chain, we get good agreement with the experimental data for 
halo cases. 
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